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ABSTRACT 

In  this  paper,  we  discuss  tlie  use  of  a  sensitivity  equation  method  to  compute  derivatives 
for  optimization  based  design  algorithms.  I’he  problem  of  de.signing  an  optima!  forebody 
simulator  is  used  to  motivate  the  algorithm  and  to  illustrate  the  basic  ideas.  Finally,  we 
indicate  how  an  existing  CFD  code  can  be  modified  to  compute  .sensitivities  and  a  nunx'rical 
example  is  presented. 
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1  Introduction 

A  large  number  of  identification,  control  and  design  problems  may  be  formulated  as  infinite 
dimensional  optimization  problems.  These  problems  arise  in  almost  all  fields  of  scieiice 
and  engineering  and  range  in  scope  from  inverse  problems  in  seismology,  to  LQK  and 
control,  to  shape  optimization  in  fluid/structure  dynamics.  See  [4,7]  for  typical  applica¬ 
tions.  Although  there  are  numerous  approaches  to  solving  these  problems,  each  approach 
requires  that  some  type  of  approximation  be  introduced  at  some  point  in  the  design  pro¬ 
cess.  Moreover,  it  is  often  the  case  that  some  iterative  scheme  is  needed  to  .solve  the  state 
equations  (in  black-box  methods  [3,8])  and  the  adjoint  equations  (in  adjoint  and  “one-shot" 
methods  [9]).  Also,  the  optimization  algorithm  may  itself  be  iterative.  In  any  case,  the 
development  of  computational  methods  for  optimal  design  and  control  can  produce  several 
levels  of  approximation  and  the  convergence  properties  of  the  overall  algorithm  are  very 
much  dependent  on  the  approximations.  In  this  paper  we  concentrate  on  the  problem  of 
computing  accurate  .sensitivities  for  gradient  based  optimization  algorithms.  In  order  to 
keep  the  paper  short  and,  at  the  same  time  illustrate  the  basic  idea,  we  concentrate  on  a 
particular  application.  We  give  a  brief  description  of  the  problem  and  use  this  problem  to 
motivate  the  algorithm  presented  below. 


2  Optimal  Design  of  a  Forebody  Simulator 

This  problem  is  a  2D  version  of  the  problem  described  in  [1,2,6].  The  Arnold  Engineering 
Development  Center  (AEDC)  is  developing  a  free-jet  test  facility  for  full-scale  testing  of 
engines  in  various  free  flight  conditions.  Although  the  test  cells  are  large  enough  to  house 
the  jet  engines,  they  are  too  small  to  contain  the  full  airplane  forebody  and  engine.  Thus, 
the  effect  of  the  forward  fuselage  on  the  engine  inlet  flow  conditions  must  be  “simulated.'’ 
One  approach  to  solving  this  problem  is  to  replace  the  actual  forebody  by  a  smaller  object, 
called  a  “forebody  simulator”  (FBS),  and  determine  the  shape  of  the  FBS  that  produces 
the  best  flow  match  at  the  engine  inlet.  The  2D  version  of  this  problem  is  illustrated  in 
Figure  1  {see  [1,2,5, 6]). 

The  underlying  mathematical  mode!  is  based  on  conservation  laws  for  mass,  momentum 
and  energy.  For  inviscid  flow,  we  have  that 
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The  velocity  components  u  and  v,  the  pressure  P,  the  temperature  T,  and  the  Mach  number 
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At  the  inflow  boundary,  we  want  to  siinnlale  a  free-pt.  so  ihai  we  sjjeciry  ttie  loial 
pressure  /b.  the  total  teinperat lire  7ii.  and  the  Maeli  nunilx'r  A/i.  W'e  also  set  !  =:  0  at 
the  inflow  boundary.  If  u/.  Pi,  and  Ti  denote  the  inflow  values  nf  tiie  r-coinponeui  of  t  lie 
velofity.  the  pressure.  an<l  tlie  teniperatnre.  llie.se  jiiay  bf'  rerovered  from  'I',,,  p.  and  A/  bv 
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The  components  of  Q  at  the  inflow  may  then  be  determined  from  (  t)  thrf)up;h  th<'  relations 
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The  forebody  is  a  solid  surface,  so  that  the  normal  component  of  tlie  velocity  vanishes. 

Of., 

nn\  +  vnn  —  0  on  the  forebody  .  (lij 

where  n\  and  112  are  the  components  of  the  unit  normal  vector  to  the  boundary.  Note  that 
we  impose  (6)  on  the  velocity  components  u  and  r.  and  not  on  the  momentum  components 
m  and  n.  Insofar  as  the  state  is  concerned,  it  is  clear  that  it  does  not  make  any  difference 
whether  (6)  i.s  imposed  on  rn  and  n  or  on  u  and  c,  since  m  =  pu  and  u  =  pi-  and  p  ^  0.  It 
can  be  shown  that  it  does  not  make  any  difference  to  the  sensitivities  a.s  well. 

.Assume  that  a,t  x  =  3  the  desired  steady  state  How  Q  =  Q{i/)  is  given  as  data  on  the 
line  (called  the  Inlet  Reference  Plane} 


IR  ^  ~  {{x,  i/)|x  =:  /f,  e  <  y  <  «^}. 

Also,  we  assume  here  that  the  inflow  (total)  Mach  number  Wo  can  b"  used  a<;  a  design 
(control)  variable  along  with  the  shape  of  the  forebody.  Let  the  forebody  be  determined 
by  the  curve  T  =  r(x),«  <  x  <  /?  and  let  p  =  (Mo,r(  )).  The  problem  ran  be  stated  a.s  the 
following  optimization  problem; 

Problem  FBS.  Given  data  Q  =  Q(i/)  on  the  IRP,  find  the  parameters  p*  =  r*(  )) 

such  that  the  functional 

J{p)  ^  \  j'  \\Qc.{0,y)  -  Q(yWdy 


2 


is  niinimized,  where  Qrxt'''.,v)  ==  solution  to  tlie  steadv  stale  Kuler  iMpiaixui 

G(Q.p)=  t^F,  +  ^F,.  =0. 

ih-  (hj 

Clearly  the  statement  oftlie  problem  is  not  eoinplele.  I'or  example,  one  slmuid  i  arefnllv 
specify  the  set  of  admissible  curves  !’(•)  and  ipiestions  remain  about  existetice.  uniquem'S' 
and  inte^raliilily  of  “the”  solution  Qx.-  We  will  not  address  tiies(‘  issues  in  this  short  note. 

Most  o])tinuzation  based  dc'sign  methods  re<piire  the  computat  ui  of  the  derivative> 
TiyQxdJ'.  (/.  ?>)■  These  derivaliv(>s  are  called  sensitivitie.s  and  various  scheme's  have  been  de¬ 
veloped  to  aj>proxjmate  tlie  sensitivities  numerically  (s('e  [d.o.lO,!  1]).  A  common  ap])roa(  h 
is  to  use  finite  differences.  In  particular,  the  steady  slate  equation  (-S)  is  solved  for  p  am! 
airain  for  p  +  Ap  and  then  ^Qxf-r.p  p)  is  approximated  by 

method  is  often  costly  and  can  introduce  large  errors.  .Another  apjrroacli  is  to  first  derive 
an  equation  (the  sensitivity  equation)  for  Q'  =  ;/.]>)  and  tlnm  numerically  solve  thi' 

equation.  VVe  shall  ilbi.strate  this  ai)proach  for  the  forebody  design  ]>roblem  and  present  a 
comparison  of  the  two  methods. 


3  Sensitivities  with  Respect  to  the  Inflow  Mach  Number 

First,  we  consider  the  design  parameter  A/J.  Thus,  we  will  derive  equations  for  the  sensi¬ 
tivity 
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The  differential  equation  system  (i)  has  no  explicit  dependence  on  the  design  param¬ 
eter  ,'V/|,  so  that  equations  for  the  components  of  Q'  are  easily  determined  by  formally 
differentiating  (1)  with  respect  to  A/J.  The  re.siilt  is  the  .sy.sten) 
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and  where,  through  (3),  the  sensitivities  (8)  and  (11)  are  related  by 
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Note  that  (9)  is  of  the  same  form  as  ({),  with  a  ciitferent  fiiix  vector,  in  jiart i(  ular,  (9} 
is  in  coiLservation  form.  /\s  a  result  of  the  fact  that  (9)  is  Initui-  in  the  juiiiied  varial)les. 
and  that  by  (12)  u',  v'.  and  P'  are  linear  in  the  component.-,  of  Q',  (9i  is  a  linear  sv.stem  in 
the  sensitivity  (7),  t.e.,  in  the  components  of  Q'. 

Now,  we  need  to  discuss  the  boundary  conditions  for  Q'.  I’,\'(e[)t  for  the  inliow  ((jiidi- 
tions,  all  boundary  conditions  are  independent  of  the  design  paramet<'r  .1/,;.  Ihiis.  the  latter 
may  be  differentiated  with  respect  to  Mq  to  obtain  boundary  conditions  for  the  sensitivities. 
For  example,  at  the  forebody  where  (6)  holds,  we  simply  would  have  that 

it'rii  +  v'Ti2  =  0  on  the  forebody  .  ( Kl) 


Similar  operations  yield  boundary  conditions  for  the  sensitivities  along  symmetry  lines, 
other  solid  surfaces,  and  at  the  outflow  boundary.  Note  that  if  instead  of  (6),  one  interprets 
the  no  penetration  condition  as  one  on  the  momentum,  i.e.,  mui  +  nn->  -  0  on  the  forebody, 
then  instead  of  (13)  we  would  have  that 

m'ni  +  n'7i2  =  0  on  the  forebody  (M) 

which  is  seemingly  different  from  (13).  However,  (6)  and  (12)  can  be  used  to  show  that 


m'm  +  n‘n2  —  p{u'ni  -f  v'lii)  +  p'{nn\  4-  vus)  =  p(u'n)  +  I'v.^ ) 


(lo) 


so  that,  since  p  ^  0,  (13)  and  (14)  are  identical. 

The  inflow  boundary  conditions  for  the  sensitivities  may  be  determined  by  differentiat¬ 
ing  (4)  and  (5)  with  respect  to  the  design  parameter  Note  that  thi.s  paratneter  appears 
explicitly  in  the  right-hand-sides  of  the  equations  in  (4)  and  (b).  Without  dilfictiity,  one 
finds  from  (15)  that 
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4  Sensitivities  with  Respect  to  the  Forebody  Design  Parameters 

We  as.sume  that  the  forebody  is  described  in  terms  of  a.  finite  number  of  dicsign  jiaraniefer?, 
which  we  denote  by  k  =  1, . . A',  and  that  the  forehody  may  be  descrilx'd  by  the  r<'latioii 

y  =  <P(x\  f'l ,  IN,  •  •  ■ ,  /'k  ).  e  <■‘•<‘'1  ( 1  ^ ) 

We  express  the  dependence  of  the  state  varialile  Q  on  the  coordinates  and  the  ilesign 

parameters  by  Q  =  Q{t,r,y,  M^,  Pi,  P'2 _ _  Pk)-  VV'e  liave  already  s('en  what  e(|nalions  can 

1)('  used  to  det(>rmine  the  sen.silivity  of  the  state  with  resix'cl  to  M~,  i.i..  for  Q'.  We 


now  discuss  what  equations  can  be  used  to  ilctenniiie  tin'  sensitivities  with  respect  tu  the 


forebod  V 


parameters  k  =  1, . . , ,  A',  i.e.,  for 
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and  E'k  = 


k  =  \ . K.  (20) 


System  <  l)  has  no  explicit  dependence  ru  the  design  parameters  Z'*.  so  that  ecjuations 
for  the  components  of  Qjt  are  easily  determined  by  differentiating  (1)  with  respect  to  P*. 
k  =  1 _ _  A'.  This  produces  the  systems,  i  =  1, . . . ,  A',  given  by 

dQk  dFk,  dFk-j 
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Moreover,  by  (3),  the  sensitivities  (20)  and  (23)  are  related  by 
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for  A:  =  1, , . A. 

All  boundary  conditions  except  the  one  on  the  forebody  also  do  not  depend  on  the 
fcrrebody  design  parameters  P*,  k  =  1,...,A'.  For  example,  consider  the  inflow  boundary 
conditions  (4)-(5).  Differentiating  these  with  respect  to  P*,  A-  =  1, _ K  yields  that 

Pki  —  aiki  =  iiki  =  Ek!  =  Tki  —  Pk!  —  Tiki  =  Tiki  =  0 ,  (25) 


at  the  inflow  boundary.  Now,  consider  the  boundary  condition  (6)  on  the  forebody.  We 
have  that  on  the  forebody 


(Combining  (6)  and  (26)  we  have  that 
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along  the  forehocly  or,  di.splaying  the  full  fiuiclioiial  (lepeiuieiice  on  tlie  co(inriiiat<:i 
design  parameters,  we  have  at  a  point  (r.y)  on  the  forebody,  and  at  any  time  t. 


ut^t,x,y  =  . Pk)\  ■  ■:PKj  T-(^:  P^■  P'^ . Pk) 

-  x:  {l.x,y^  <l>{a-,  P^.  ,  Pk).  MJ,  /', .  P.. .  . . .  Pk 


We  can  proceed  to  differentiate  (28)  with  respect  any  of  the  forehody  design  j)araineters 
Pk,  k  =  I, . . . ,  K .  The  result  is  that,  along  the  forebody  for  A:  -  l . A'. 

/ du\  /  \  f  d<^\  fdr\f(n\ 

dx  vdj/j  j  vdx  j  “dr  \dPk)  ^  ydvj  Vd/\  j 


fd<P\  d 
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dr  \dPk 
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where  u,  v,  and  their  derivatives  are  evaluated  at  the  forebody  (r,y  'hfr)). 

If  an  iterative  scheme  is  used  to  find  a  steady  state  solution  of  this  system  ((21). 
(25),  (29)),  then  we  as.sume  that  present  guesses  for  the  state  variable.s  u  and  c  and  their 
derivatives  du/dy  and  dv/dy  and  for  the  design  parameters  A/|  and  Pk,  k  =  1, , ,  ,  A.  are 
known.  It  follow's  that  the  right-hand-side  of  (29)  is  known  as  w-ell  and  equation  (29).  the 
boundary  conditions  along  the  forebody  for  the  sensitivities  with  respect  to  the  forebody 
design  parameters,  is  merely  an  inhomogeneous  version  of  (27).  the  boundary  condition 
along  the  forebody  for  the  state. 

Let  us  now  specialize  to  the  type  of  forebodies  considered  by  Huddleston.  [5,0].  i.e. 

K 

4>(x;P,,P2....,Pk)  =  ^Pk<t>k{x),  (50) 

k=l 

where  4)k{x),  it  =  1,...,  A',  are  prescribed  functions,  e.g.,  Bezier  curves.  In  this  case, 


=  <t>kix)  and 


\dPkJ 


Combining  (29)-(32),  one  obtains  that,  at  any  point  (z,^(a:))  on  the  forebody  and  for  each 
k  t=  I, ....  K, 


;=1 


Uk  -  Vk  - 


4>k  -  x‘-} — I- 
dx 


For  forebodies  of  the  type  (30),  (33)  gives  the  the  boundary  conditions  along  the  forebody 
for  the  sensitivities  with  respect  to  the  forebody  design  parameters  Pk,  k  ~  1,.  ,  .,A  .  It  is 
now  clear  that,  given  guesses  for  the  state  variables  u  and  v  and  their  derivatives  Pn/Py  and 

dx)/dy  and  for  the  design  parameters  Aig  and  P*,  k  ~  I . A,  then  the  right-hand-side  of 

(33)  is  known. 


5  Computing  Sensitivities  using  an  Existing  Code  for  the  State 

Suppose  one  has  available  a  code  to  compute  the  state  variables,  t.r.,  to  find  approximate 
solutions  of  ( 1 )  along  with  boundary  and  initial  conditions.  In  principle,  it  is  an  easy  matter 
to  amend  stich  a  code  so  that  it  can  also  compute  .sensitivities. 
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First,  let  iis  compare  (1)  with  (9).  If  one  wishes  to  amend  the  existinj'  state  ciide  that 
can  handle  ( 1 )  so  that  it  can  treat  (9)  as  well,  one  has  to  change  the  definitit  ns  of  the  finx 
functions  from  those  given  in  (2)  to  those  given  in  ( 10).  Note  that  the  solutioti  for  tiie  state 
is  needed  in  order  to  evaluate  the  flux  functions  of  (10). 

Next,  note  that  (9)  and  (22)  are  identical  differential  equations.  Thus,  the  changes 
made  to  the  code  in  order  to  treat  (9)  can  al.so  be  used  to  treat  (22).  In  fact,  a.s  long  a.s 
the  differentia!  equation  and  any  other  part  of  the  problem  sjtecification  do  not  explicitly 
depend  on  the  design  parameters,  the  analogous  relations  w'ill  be  the  same  for  all  the 
sensitivities. 

The  only  changes  that  vary  from  one  sensitivity  calcnlation  to  another  are  tho.se  that 
arise  from  conditions  in  which  the  design  parameters  appear  explicitly.  In  our  example,  for 
the  sensitivity  w'ith  respect  to  Mq,  one  must  change  the  portion  of  the  code  that  treats  tlie 
inflow'  conditions  (4)-(,5)  so  that  it  can  instead  treat  (16)-(17),  The  only  changes  needed 
to  accomplish  this  are  to  the  data  of  the  inflow  conditions.  In  the  prol)!em  considered 
here,  the  nature  (cc.  what  variables  are  specified)  of  the  l)oundary  conditions  at  the  inflow, 
and  everywhere  else,  is  not  affected.  Note  that  for  the  sensitivity  with  respect  to  the 
boundary  condition  (13)  on  the  forebody  is  the  same  as  that  for  the  state,  given  by  (6). 

For  the  sensitivities  with  respect  to  the  forebody  design  parameter.s,  the  inflow  boundary 
conditions  simplify  to  (25),  t.e.,  they  become  homogeneous.  The  boundary  condition  at  the 
forebody  is  now  given  by  (29)  oi  (33).  Once  again,  the  nature  of  the  boundary  conditions 
i.s  unchanged  from  that  for  the  state,  and  only  the  data  that  is  specified  is  different.  For  the 
inflow  boundary  conditions,  we  may  still  specify  the  same  conditions  for  the  sensitivities, 
but  now  they  would  be  homogeneous.  The  boundary  conditions  along  the  forebody  change 
only  in  that  they  become  inhomogeneous,  (compare  (27)  and  (33)). 

in  summary,  to  change  a  code  for  the  state  so  that  it  also  handles  the  sensitivities,  one 
must  redefine  the  flux  functions  in  the  differential  equatiors,  and  the  data  in  the  boundary 
conditions.  The  changes  necessary  in  the  code  to  account  for  any  particular  relation  that 
does  not  explicitly  involve  the  design  parameters  are  independent  of  which  sensitivity  one 
is  presently  considering. 

The  previous  remarks  are  concerned  only  with  the  changes  one  must  effect  in  a  state 
code  in  order  to  handle  the  fact  that  one  is  discretizing  a  different  problem  when  one 
considers  the  sensitivities.  We  have  seen  that  these  changes  are  not  major  in  nature. 
However,  there  are  additional  changes  that  may  be  needed  when  one  attempts  to  solve 
the  discrete  equations.  In  the  numerical  results  presented  below  we  use  the  finite  difference 
code  “PARC”  (see  [2,-5])  to  solve  the  state  and  .sensitivity  equations.  However,  the  following 
comment.^  apply  equally  well  to  other  CFI)  codes  of  this  type. 

Since  we  are  interested  in  the  steady  design  problems,  the  time  derivative  in  (!)  i.s 
considered  only  to  provide  a  means  for  marching  to  a  steady  state.  Now,  suppose  that  at 
any  stage  of  a  flanss-Newton,  or  other  iteration,  we  have  used  F’.XRC  to  find  an  approximate 
steady  stale  solution  of  (1)  plus  boundary  conditions.  In  order  to  do  this,  one  has  to  solve 
a  .sequence  of  linear  algebraic  systems  of  the  type 

(/  +  AtA(Qj,"’))Q<,’‘  +  '^=  (Qi”UA/B(Q',"')),  u  =  0,  1,2,...  ,  (.34) 

where  thf‘  sequence  is  teniiinated  when  one  is  satisfied  that  a  steady  state  has  been  reached 
and  where  q)”’  denote.s  the  discrete  approximation  to  the  slate  Q  at  the  lime  t  =  uA/.  We 
denote  this  .steady  state  .solution  for  the  apf>roximation  to  the  .state  by  Q/, .  One  prol»h'm  of 
the  ty|)e  f34)  is  solved  for  every  time  step.  In  (31).  the  matrix  A  aiul  vc'ctor  D  ari.se  from 
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thp  spatial  discretization  of  the  fluxes  and  the  houndary  condilions.  Hoili  if  the>e  dein  tui 
on  the  state  at  the  previous  time  level. 

Having  computed  a  steady  state  solution  by  (dl),  the  task  at  hand  is  to  now  ciiiiipute  i  he 
sensitivities.  We  will  focus  on  Q',  the  sensitivity  with  respect  to  the  irifiou  .Marli  nutiiher. 
Analogous  results  hold  for  the  sensitivities  with  respect  to  tiie  forehody  design  parameters. 
Recall  that  given  a  state,  the  sensitivity  equations  are  linear  in  the  sensitivities,  rherefore. 
if  one  is  interested  in  the  steady  state  sensitivities,  instead  of  (9)  one  may  directly  treat  its 
stationary  version 


(>F',  rJF'o 

^  +  ^  =  0. 

ax  ay 


f ;(.')) 


Since  (35)  is  linear  in  the  components  of  Q',  one  does  not  need  to  consider  marcliiiig 
algorithms  in  order  to  compute  a  steady  sensitivity.  One  merely  discretizes  (3.5)  ami  so!ve^ 
the  resultant  linear  system,  which  has  the  form 


A'(Qfc)Q',  =B'(Qa), 


(36) 


where  Q),  denotes  the  discrete  approximation  to  the  steady  sensitivity.  The  matrix  A'  and 
vector  B'  differ  from  the  A  and  B  of  (34)  because  we  have  discretized  different  differential 
equations  and  boundary  conditions.  Note  that  A’  and  B'  in  (36)  depend  only  on  the  steady 
state  Q/i  and  thus  (36)  is  a  linear  system  of  atgebrci,  ^quattons  for  the  discrete  sensitivity  Q)^. 

The  cost  of  finding  a  solution  of  (36)  is  similar  to  that  for  finding  the  solution  of  (34) 
for  a  single  value  of  n,  j.e.,  for  a  single  time  step.  The  differences  in  the  assembly  of  the 
coefficient  matrices  and  right-hand-sides  of  (34)  and  (36)  are  minor.  Thus,  in  theory  at 
least,  one  can  obtain  a  steady  sensitivity  in  the  same  computer  lime  it  lakes  to  perform  one  time  step 
in  a  state  calculation.  If  one  wants  to  obtain  all  the  sensitivities,  e.g..  A' +  i  in  our  example, 
one  can  do  so  at  a  cost  similar  to,  e.g.,  K  1  time  steps  of  the  state  calculation.  This 
is  very  cheap  compared  to  the  multiple  state  calculations  necessary  in  order  to  compute- 
sensitivities  through  the  use  of  difference  quotients. 

In  practice,  these  “optimal”  estimates  of  speed  up  are  rarely  achieved.  Moreover,  it  is 
important  to  note  that  finite  difference  (FD)  and  sensitivity  equation  (SE)  methods  do  not 
necessarily  produce  the  same  re.sults.  Since  the  ultimate  goal  is  to  find  useful  and  cheap 
gradients  for  optimization,  the  most  important  is.sue  is  whether  or  not  the  SE  method 
combined  with  an  optimization  algorithm  produces  a  convergent  optimal  design  as  fast  as 
possible.  We  have  tested  this  scheme  on  the  forebody  design  problem  with  excellent  results. 


6  A  Numerical  Example 

In  order  to  illustrate  the  use  of  the  SE  method  in  computing  sensitivities,  we  used  the 
PARC  code  as  described  above  to  find  appro.ximate  solutions  of  the  sensitivity  equations 
and  compared  the  results  to  the  finite  difference  method.  In  Figure  2  we  show  the  ap¬ 
proximations  of  ,  Pi,  P2)  for  Mo  =  2, /’i  —  .1  and  Po  -  -15  where  the  forebody 

is  described  by  two  Bezier  parameters  (Pi,P2)-  Both  pictures  are  “converged”  estimates. 
Note  that  there  are  considerable  differences  between  the  FD  method  and  the  SE  method. 
Moreover,  in  Figure  3  we  see  that  not  only  do  the  FD  and  SE  methods  produce  different 
sensitivities,  the  value  of  the  step  size  APi  can  greatly  influence  the  FD  approximatioius. 
Finally,  we  note  that  the  SE  method  ran  4  to  .5  tiine.s  faster  than  the  FD  method.  .-\lso, 
although  '  pace  prohibits  a  discussion  of  the  optimization  problem  here,  we  have  u.sed  lh(' 
SE  method  in  a  trust  region  optimization  scheme  to  produc*'  an  optimal  forebody  design 
for  the  21)  problem  in  [5,6].  1  hose  results  will  appttar  in  a  foithronung  jiaper. 
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7  Conclusions 

rhc  problem  of  conii)«ting  acniral-,'  sensitivities  in  pr(ibl<’ins  involving  solution^,  to  paiain 
eteri?.(Hi  partial  (iilFerentiai  eqiialioiis  is  an  iiiiportaril  {)art  of  <)[)!imal  desioti.  The  f>oa!  is  to 
fiml  derivatives  of  solutions  of  partial  differential  equations  witli  respei  t  to  various  param¬ 
eters  (iiuhiding  domain  shapes)  and  to  use  these  derivatives  in  some  type  of  optimization 
scheme.  In  almost  all  practical  problems,  solutions  must  be  obtained  by  numerical  ajiprox- 
imations.  This  fact  lead.s  to  “black  bo.v”  methods  for  optimal  design.  In  its  nio.st  basic 
form,  a  black  bo.x  method  produces  appro.ximate  solutions  that  are  then  differentiated  (by 
finite  differences,  automatic  differentiation,  etc.),  T  he  sensitivity  eciualion  (SE)  method 
presented  here  is  based  on  first  deriving  partial  differential  equations  for  the  derivatives 
ami  then  approxitnating  these  equations  numerically.  Both  approaches  produce  numerical 
airproximations  of  the  sensitivities.  However,  the  (SE)  method  can  often  reduce  co’n|)uta- 
tional  effort,  speed  up  the  calculations  and,  p.ovided  that  accurate  computational  scheiiies 
can  be  devised  for  the  sensitivity  equations,  the  derivatives  can  lie  computed  witli  the  satin 
degree  of  accuracy  as  the  state.  The  2D  op  '.mal  forebody  simulator  problem  is  an  excellent 
problem  for  illustrating  these  points.  The  numerical  results  presented  here  show  that  the 
(SE)  method  is  potentially  applicable  to  real  problems  and,  at  the  same  time,  raises  many 
interesting  theoretical  and  practical  questions. 
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